The Structure of the Pairing Interaction in the 2D Hubbard Model 
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Dynamic cluster Monte Carlo calculations for the doped two-dimensional Hubbard model are used 
to study the irreducible particle-particle vertex responsible for d^2_y2 pairing in this model. This 
vertex increases with increasing momentum transfer and decreases when the energy transfer exceeds 
a scale associated with the Q — {n, n) spin susceptibility. Using an exact decomposition of this 
vertex into a fully irreducible two-fermion vertex and charge and magnetic exchange channels, the 
dominant part of the effective pairing interaction is found to come from the magnetic, spin S = 1 
exchange channel. 



Numerical studies of the doped two-dimensional Hub- 
bard model have shown that strong d^2_y2 pairing corre- 
lations develop as the temperature is lowered and re- 
cent work has provided evidence that the Hubbard model 
near half filling does have a superconducting ground state 
0, 0, Inspite of this progress, the goal of using nu- 
merical methods to determine the nature of the pairing 
mechanism has proven elusive. In this paper we de- 
velop a new approach which combines numerical and di- 
agrammatic methods to determine the structure of the 
pairing interaction in the doped two-dimensional Hub- 
bard model. Our study will focus on the 4-point ver- 
tex calculated with a Quantum Monte Carlo dynamic 
cluster approximation 5, 6, JJ. From this vertex and the 
QMC/DCA results for the single particle Green's func- 
tion, we have determined the irreducible particle-particle 
and particle-hole vertices. The leading low tempera- 
ture eigenvalue of the Bethe-Salpeter equation for the 
particle-particle channel is shown to have dx2_y2 sym- 
metry. We then examine the momentum and energy de- 
pendence of the irreducible particle-particle vertex. De- 
composing this vertex into the sum of a fully irreducible 
two-fcrmion vertex and particle-hole exchange magnetic 
[S = 1) and charge density (S* = 0) channels, we find 
that the dominant contribution to the pairing interac- 
tion comes from the magnetic (5 — 1) exchange. 

The Hubbard model that we study has a near- 
neighbor, one-electron hopping t and an on-site Coulomb 
interaction U . 

H = -tY^ {c\^Cj„ -f c]^c,^) + U ^n^tna 

We will take U/t — 4 and adjust ^ so that the average site 
occupancy (rii) = 0.85. We have carried out dynamical 
cluster Monte Carlo calculations IH IE for the 24-site 
/c-cluster shown in the inset of Fig. 2. 



For a 2D system the dynamical cluster approximation 
maps the original lattice model onto a periodic cluster of 
size Nc — embedded in a self-consistent host 0, S ■ 
The essential assumption is that short-range quantities, 
such as the self energy and its functional derivatives (the 
irreducible vertex functions) are well represented as dia- 
grams constructed from the coarse-grained Green's func- 
tion. To this end, the first Brillouin zone is divided into 
Nc cells, with each cell represented by its center wave- 
vector K surrounded by N/Nc lattice wavevectors labeled 
by k. The reduction of the A^-site lattice problem to an 
effective Nc site cluster problem is achieved by coarse- 
graining the single-particle Green's function, i.e. averag- 
ing G(K + k) over the k within a cell which converges 
to a cluster Green's function Gc(K). Consequently, the 
compact Feynman diagrams constructed from Gc(K) col- 
lapse onto those of an effective cluster problem embedded 
in a host which accounts for the fluctuations arising from 
the hopping of electrons between the cluster and the rest 
of the system. The compact cluster quantities are then 
used to calculate the corresponding lattice quantities. 

For example, the DCA cluster one- and two-particle 
Green's functions that we calculate have the standard 
finite temperature definitions 



(2a) 



and 



Gc2(y4.---cn{X4:, X:^; X2, Xi) — 

- {T,c,jX,)c,jX,)cl{X2)cl{X,)) . (2b) 



Here, Xi — pi.i,Te), where denotes a site in the 
DCA cluster, re the imaginary time, Tr is the usual r- 
ordering operator, and Ca \X2) creates (destroys) a par- 
ticle on the cluster with spin a. Fourier transforming on 
both the cluster space and imaginary time variables gives 
Gc{K) and Gc2{Ki, Ky, K2, Ki) with K = {K,iLJn,a). 
Using Gc{K) and Gc2{Ki,K^]K2,Ki), one can extract 



2 





r 






+ 






r 








r 






+ 


r"" 




r 








r"" 






+ 




+ 












momentum fc' and Matsubara a;„' variables. We decom- 
pose k' = K' + k'. By assumption, irreducible quanti- 
ties like TPP and (j)a do not depend on k', allowing us to 
coarse-grain the Green function legs, yielding an equation 
that depends only on coarse-grained and cluster quanti- 
ties 



FIG. 1: Bethe-Salpeter equations for (a) the particle-particle 
and (b) the particle-hole channels showing the relationship 
between the full vertex, the particle-particle irreducible ver- 
tex r'^'', and the particle-hole irreducible vertex F^^, respec- 
tively, (c) Decomposition of the irreducible particle-particle 
vertex F^^ into a fully irreducible two-fermion vertex Airr plus 
contributions from the particle-hole channels. All diagrams 
represent DCA cluster quantities, including the Green func- 
tion legs. 



the cluster four-point vertex F from 
Gc2{Ki, K3; K2, Ki) = 

— Gc{Ki)Gc{K2) [5Ki,Ki&K2,K3 ~ ^Ki,K3&K2,K4\ 

T 

SKi+K2,K3+KiGc{K4) Gc{K3)r {K4, K3; K2, Ki) 



N 

X Gc{K2)Gc{Ki). 

Then, using Gc and F, one can determine the irreducible 
particle-particle and particle-hole vertices Fpp and FP^ 
from the Bethe-Salpeter equations shown in Fig.QIa) and 
(b) . There is a second particle- hole vertex but it is simply 
related to T^^. Note that FP*^ and Fpp don't have a c 
suffix, since both the lattice and the cluster share these 
compact quantities. Because of the rotational invariance 
of the Hubbard model, it is convenient to separate the 
particle-particle channels into singlet and triplet and the 
particle-hole channels into a magnetic part which carries 
spin S — 1 and a charge density part which has S* = 0. 

In order to determine the nature of the low tempera- 
ture correlations, we use these irreducible vertices and the 
lattice single particle Green's function to calculate the 
Bethe-Salpeter eigenvalues and eigenfunctions in various 
channels. For example, in the particle-particle channel 
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^FPP iK,-K;K',-K') 



\MK) (5) 



with ^^\K') = f Ek' Gt(K' + k') Gi(-K' - k') 
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X GT(fc')Gi(-fc')0a(i^') = Aa0a(if) (4) '^d 



with a similar equation using FP*^ for the particle-hole 
channel. Here, the sum over k' denotes a sum over both 



T/t 



FIG. 2: Leading eigenvalues of the Bethe-Salpeter equation 
(e.g. Eq. 13 in various channels for U/t = A and a site occu- 
pation {n) — 0.85. The Q = (tTjTt), LOm = 0, 5 = 1 mag- 
netic eigenvalue is seen to saturate at low temperatures. The 
leading eigenvalue in the singlet Q = (0, 0), Um ~ particle- 
particle channel has d^2_y2 symmetry and increases toward 
1 at low temperatures 0- The largest charge density eigen- 
value occurs in the Q = (0, 0), iOm = channel and saturates 
at a small value. The inset shows the distribution of fc-points 
for the 24-site cluster we have studied. 



In Fig. 121 we show the leading eigenvalue versus tem- 
perature for the pairing, charge density, and magnetic 
channels for U/t = 4 and (n) = 0.85. As the temperature 
is reduced, the leading particle-hole eigenvalue occurs in 
the magnetic channel and has a center of mass momen- 
tum Q = (tt, tt) and ojm = 0. Previous Monte Garlo 
calculations on 8x8 lattices show that for this doping 
the peak response is slightly shifted from (tTjTt), but our 
24-site cluster lacks the resolution to see this This 
antiferromagnetic eigenvalue grows and then saturates 
at low temperatures. The leading particle-particle eigen- 
value is a spin singlet and, as shown in the inset of Fig. 3, 
its eigenfunction /pd 2_ 2 dx^-y^ symmetry. The tUn 
frequency dependence of the normalized gap function 
(K, LUn) at the antinodal point K — (tt, 0) is plot- 
ted in Fig. 3. As shown, it is even in w„, corresponding to 
a dx2-y2-weLve singlet, even frequency pairing. Also plot- 
ted in this figure is the Wm-dependence of the Q = (tt, tt) 
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spin susceptibility x(Q, Wm) normalized to coincide with 
4>d^2_ 2(K,t<;„) at (jJn = ttT. The Matsubara frequency 
which enters the gap function corresponds to a fermion 
frequency ijJ„ — (2n + l)7rT', while = 2mnT for the 
spin susceptibility, leading to the interlacing of points 
shown in Fig. 3. From the momentum and frequency de- 
pendence of the gap function (/)(K, ujn), it follows that the 
irreducible particle-particle vertex is an increasing func- 
tion of the momentum transfer and is characterized by 
the same energy scale that enters the spin susceptibility 
x(Q,Wm)- At larger values of U it will be interesting to 
see if there is an increased tendency for a finite response 
at large Matsubara frequencies indicating a contribution 
from the upper Hubbard band. 




FIG. 3: The Matsubara frequency dependence of the eigen- 
function (pd 2 2 (-^i leading particle-particle eigen- 

value of Fig. 2 for K = (tt, 0) normalized to (p(K, ttT) (red, 
solid). Here, = (2n + 1)-kT with T = 0.125f. The Mat- 
subara frequency dependence of the normalized magnetic spin 
susceptibility 2x(Q, t^™)/[x(Q, 0)+x(Q, ^ttT)] for Q = (tt, tt) 
versus = 2m7rT (green, dashed). Inset: The momen- 
tum dependence of the eigenfunction (j)d 2_ 2 (K, ttT) normal- 
ized to 2 2 ((0, 1"), 7i"r) shows its d^2_y2 symmetry. Here, 
oJn = ttT and the momentum values correspond to values of 
K which lay along the dashed line shown in the inset of Fig. 2. 



To learn more about the mechanism responsible for 
dx2-y2 pairing in the doped Hubbard model, it is use- 
ful to decompose the pairing interaction PPP as shown in 
Fig. Ic. Here, the irreducible particle-particle vertex is 
given as a combination of a fully irreducible two-fermion 
vertex Airr andpartially reducible particle-hole exchange 
contributions 0, llOj . For the even frequency, even mo- 
mentum part of the irreducible particle-particle vertex 
TPP^^{K,K') = l/2[TPPiK,K') + TPPiK,-K% we ob- 
tain 



with K = (K, iujn] 
charge density {S 
hole channels 

1 



The subscripts d and m denote the 
= 0) and magnetic (S* = 1) particle- 



(7) 



T^/.^{K - K'- K', -K) - T^I]^{K - K'; K', -K) 



Ta,^m{K + K'- -K', -K) - T^^JK + K'; -K', -K) 



d/n 
d/n 



-MK,K') + -<i>,r,{K,K') (6) 



The center of mass and relative wave vectors and frequen- 
cies in these channels are labeled by the first, second and 
third arguments respectively. Using the Monte Carlo re- 
sults for G and F, we have solved the matrix equations 
shown in Figs. 1(a) and (b) to determine Fpp, and 
$,„. Then, substituting these into Eq. ©, we have de- 
termined the fully irreducible vertex Airr. 

Monte Carlo results for the irreducible particle-particle 
vertex Fpp obtained from the 24-site cluster approxima- 
tion are shown in Fig. 4a. Here, we set w„ = lo'^ — vrT, 
K = (tt, 0) and K' takes momentum values along the 
dashed line shown in the inset of Fig. 2. As the tempera- 
ture is lowered, Fpp increases as the momentum transfer 
q = K — K' increases as one expects for a d-wave pair- 
ing interaction. To understand the origin of this behav- 
ior, the contributions of the particle-hole {S = 0) charge 
density and (5* = 1) magnetic channels are plotted in 
Figs. 4c and 4d respectively and the contribution from 
the fully irreducible vertex Airr is shown in Fig. 4b. It is 
clear that the dominant contribution to Fpp comes from 
the 5 = 1 magnetic channel. The charge density channel 
and the fully irreducible vertex are basically flat in mo- 
mentum and change relatively little as the temperature is 
reduced. Thus, based upon the decomposition of the ir- 
reducible particle-particle interaction shown in Fig. 4, we 
conclude that the pairing mechanism in the doped two- 
dimensional Hubbard model is mediated by the exchange 
of 5 = 1 particle- hole spin- fluctuations. 

To summarize, we have studied the pairing interaction 
FPP of a doped (n) — 0.85, two-dimensional Hubbard 
model with U/t = 4. We found that the eigenfunction 
(/)(K, iuJn) of the leading low temperature eigenvalue in 
the particle-particle pairing channel is an even frequency 
singlet with dx2_y2 symmetry. The momentum and fre- 
quency dependence of 0(K, iuJn) imply that Fpp increases 
as the momentum transfer q = K — K' increases and 
that its dynamics is set by the same characteristic en- 
ergy scale as the spin susceptibility. It was also found 
to increase as the temperature was lowered, saturating 
when the leading antiferromagnetic eigenvalue stopped 
growing. Finally, using an exact decomposition of Fpp, 
we showed that the dominant contribution to this inter- 
action comes from the S — \ particle-hole channel. We 
believe that the calculation and analysis of the four-point 
vertex provides a useful, unbiased method for determin- 
ing the nature of the leading correlations of interacting 
many-electron systems and the structure of the mecha- 
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FIG. 4: (a) The irreducible particle-particle vertex T^''' versus q = K — K' for various temperatures with Un = lu„i — nT. Here, 
K — (tt, 0) and K' moves along the momentum values of the 24-site cluster which lay on the dashed line shown in the inset of 
Fig. 2. Note that the interaction increases with the momentum transfer as expected for a d-wave pairing interaction, (b) The 
q-dependence of the fully irreducible two-fermion vertex Airr- (c) The q-dependence of the charge density (5* = 0) channel |4>d 
for the same set of temperatures, (d) The q-dependence of the magnetic {S = 1) channel |$m- 



nisms responsible for them. 
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